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Abstract 

The diffusion radius of the IS muonic hydrogen atoms in gaseous H2 targets with various deu- 
terium admixtures has been determined for temperatures T = 30 and 300 K. The Monte Carlo 
calculations have been performed using the partial differential cross sections for p/i and da atom 
scattering from the molecules H2, HD and D2. These cross sections include hyperfme transitions 
in the muonic atoms, the muon exchange between the nuclei p and d, and rotational-vibrational 
transitions in the target molecules. The Monte Carlo results have been used for preparing the 
time-projection chamber for the high-precision measurement of the nuclear u~ capture in the 
ground-state pu atom, which is now underway at the Paul Scherrer Institute. 
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Theoretical studies of the muonic atom diffusion in molecular hydrogen-isotope targets 
are important for many experiments in low-energy muon physics. In particular, knowledge 
of the diffusion radius of muonic hydrogen atoms in gaseous H-D targets is required for 
investigations of the fi~ nuclear capture in the pfi and dfi atoms created in H-D targets. The 
diffusion radius i?diff is defined as the distance between the point of the muon stop in H-D 
and the point of the muonic atom disappearance due to the muon decay or to the muon 
nuclear capture. Since the fi~~ capture rate on p or d is several orders of magnitude lower than 
the muon decay rate, Rdm is practically determined by the point of the muon decay. A high- 
precision measurement of the rate A s for the muon capture pfi — > + n in the jrround-state 
pfi atom (MuCap experiment) is underway at the Paul Scherrer Institute [H 0, H 0] . The 
rate A s for the singlet state F = of the total muonic atom spin F is sensitive to the weak 
form factors of the nucleon, especially to the induced pseudoscalar coupling constant g P . 
As a result, this experiment will provide a rigorous test of theoretical predictions based 
on the Standard Model and low-energy effective theories of QCD 0, 0| . A high-precision 
measurement of the fi~ capture rate in the process dp, — > + n + n is under consideration 
by the MuCap collaboration Such an experiment would be uniquely suited to study the 
axial meson exchange currents in the two-nucleon system. 

In this paper, main results of the Monte Carlo simulations for determining the optimal 
conditions for the MuCap experiment are presented. The time-projection chamber is filled 
with almost pure H2 gas which, however, contains a very small D2 contamination. In the 
isotope exchange process pfi + d — > dfi + p, the energy of about 135 eV is released in the 
centre-of-mass system. Therefore, the created dfi atom gains the collision energy of a few 
tenseV. As a result, the diffusion radius is significantly enlarged. This leads to an en- 
hanced absorption of the muons in the time-projection-chamber walls and limits the spatial 
resolution. The determination of the highest acceptable D2 contamination has been one of 
the aims of the presented simulations. Since the capture rate A s depends strongly on the 
total pfi spin, it is necessary to calculate the time evolution of the population of the pfi 
spin states. The initial distribution of the spin states F — 1 and F = is statistical. The 
simulations have been performed for the target temperatures T = 30 and 300 K. The target 
density has been fixed at the constant value <fi = 0.01 (relative to the liquid hydrogen den- 
sity of 4.25 x 10 22 atoms/cm 3 ), which corresponds to the pressure of about 9 bar at 300 K. 
At such a density, the probability of formation of the muonic molecule ppfi is small. In 
higher-density targets, the muon nuclear capture inside ppfi is significant. This leads to 
serious problems with interpreting the experimental data owing to inaccuracy of the rate 
for the ortho-para conversion of the ppfi molecules 0, 0, 0, 0, m ■ The spin-flip transition 
Pfi(F = 1) + p — > pfi(F = 0) + p due to the muon exchange between the protons is still 
sufficiently strong at ~ 0.01 to ensure a fast quenching of the higher hyperfine state F = 1 
and, therefore, an unambiguous A s measurement. 

The Monte Carlo kinetics code includes the muon decay, pfi and dfi scattering from 
the molecules H 2 , HD and D 2 , and formation of the molecules ppfi, pdfi and ddfi. In the 
scattering process, the atoms can change their spin states. The isotope exchange reaction 
pfi + d — > dfi + p in pfi scattering from HD and D 2 is taken into account. Also, all possible 
rotational and vibrational transitions in the target molecules are included. At the collision 
energies e < 10 eV (in the laboratory system), the scattering processes are described using 
the differential cross sections da/dQ for scattering from the hydrogenic molecules 0, El 
("molecular" cross sections). At higher energies, effects of the molecular binding and electron 
screening can be neglected and, therefore, the differential cross sections for the muonic atom 



2 



scattering from hydrogen-isotope nuclei are used 0, [1(1 U , 13 ( "nuclear" cross sections) . In 
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FIG. 1: Transport (dashed lines) and total (solid lines) cross sections for the scattering oipfi(F = 0) 
and dfi atoms from a ground-state H2 molecule versus the collision energy e in the laboratory sys- 
tem. The doubled total cross sections (dash-dotted lines) for the corresponding nuclear scattering 
are shown for comparison. The hyperfine-transition threshold is denoted by AE^. 



fig. ^ the total molecular cross sections for pp(F = 0) and dp scattering from the ground- 
state H 2 molecule are shown as an example. The muonic atom spin is conserved in the 
presented processes. Also, the corresponding transport cross sections, defined as 

a tam = J dfi(l-costf)^^, (1) 

are shown. The scattering angle is denoted here by The doubled total nuclear cross 
sections for the processes pp(F = 0) + p — > pp(F = 0) + p and dp, + p — > dp + p are 
plotted for comparison. The transport cross sections approach the total cross section only 
at e — > 0, which demonstrates strong anisotropy of the molecular cross sections. Large 
differences between the molecular and nuclear cross sections at e < 1 eV are due to molecular- 
binding and electron-screening effects. The total molecular and nuclear cross sections for all 
combinations of the three hydrogen isotopes are presented in ref (l3j. 

The time evolution of the hyperfine states, the energy distribution of the muonic atoms, 
and the radial distribution of the muon decays were calculated for various initial conditions. 
All the presented results are given for a fixed target density tp = 0.01. The initial distribution 
of the pp or dp kinetic energy was described by the two Maxwell components: thermal (50%) 
and energetic (50%) with the mean energy e avg = 1-5 eV, according to the experimental 
results |l4lll5|. The calculated time evolution of the F = state and of the mean pp(F = 0) 
kinetic energy are shown in fig. for a pure H2 at T = 300 K. The pp atoms starting at 
e ~ 1 eV are slowed down within a few tens ns to energies where the spin-flip transitions 
F = — > F = 1 are impossible. The hyperfine-transition threshold is = 0.182 eV in 

the pp + p centre-of-mass system. After this time, the F = 1 state disappears with a time 
constant of 6 ns. Hence, about 50 ns after the muon stop, the relative population of the 
F = 1 state is below 0.01 and the measurement is no longer distorted by the population 
of the upper hyperfine level. All that takes place when most of the initially energetic 
atoms remains epithermal (e 3> k^T, where ks is the Boltzmann constant). The pp(F = 0) 
thermalization from e ~ 0.1 eV takes about 400 ns. As it is illustrated in fig- El the pp(F = 0) 
energy spectrum is epithermal for times much longer than in the case of pp{F — 1) atoms. 
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FIG. 2: Time dependence of the p[i(F = 0) population and of the mean pfi(F = 0) kinetic 
energy e avg in a pure H2 at T = 300 K and (p = 0.01. 
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FIG. 3: Energy distribution of p[i(F = 0) and p[i(F 
several moments after the muon stop. 
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1) atoms in a H2 gas at T = 300 K, for 



Only after the total deexcitation of the F — 1 level, the p[t(F = 0) energy distribution takes 
the final Maxwellian form with e avg = 0.04 eV. Most of the p[t diffusion until the muon 
decay takes place after the system has been thermalized. 

The mean diffusion range, which is important for the optimisation of the pressure and 
temperature of H2 filling the time-projection chamber, equals about 1 mm. However, long- 
lived (t > 10 [is) muons travel much farther, which limits the reachable spatial resolution. 
Figure 0] shows the fraction of the muon decays outside the diffusion radius from the point 
of p[i formation, for the target temperature of 30 and 300 K. The thermal diffusion is 
significantly reduced at 30 K. This effect is, however, limited because of the above-mentioned 
p[i acceleration in the spin-flip process. The radial distribution of the muon decays for several 
time intervals is plotted in fig. 

The calculated values of the mean diffusion radius for a pure H2 target at = 0.01 are 
summarised in table |1J The results are given for T = 30 and 300 K. The realistic two- 
Maxwell distributions of the initial p[i energies have been used. Also, the thermalized initial 
distributions of p[i atoms with the depleted F = 1 state have been employed in order to 
investigate the thermal part of the diffusion. 

A real H 2 target always contains a certain admixture of deuterium. Figure El demonstrates 
that the maximal muonic atom diffusion radius is greatly increased when the deuterium 
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FIG. 4: Fraction of the p~ decays outside the pp diffusion radius from the point of pp formation 
in a pure H2, for times t < 20 ps and T = 30 and 300 K. 
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FIG. 5: Radial distribution of the pT decays for the time intervals 0-5, 5-10 and 10— 15/xs at 
temperature T = 300 K. 

concentration of 10 _4 -10 -3 is present in the H2 target. The long-range tail in the radial 
distribution of the muon decays is due to very energetic (~ 45 eV) dp atoms formed in the 
pp collisions with deuterons. These dp atoms can travel at large distances owing to the 
deep Ramsauer-Townsend minimum in the d\x + p cross section (see fig. [TJ. Therefore, it 
is crucial to reduce the deuterium concentration to a very low level of about 10 -10~ 6 . 
The distribution of the dp diffusion radius in a pure D 2 gas is shown in fig. for T = 30 



TABLE I: The calculated mean diffusion radius of the pp atom in pure H2 targets for various initial 
conditions and ip = 0.01. 
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mean i?diff [mm] 
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FIG. 6: Radial distribution of the [i~ decays in H2 with the deuterium concentrations = 10 
and 1CT 4 , for times t < 20 ^s and T = 300 K. 
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FIG. 7: Radial distribution of the /i decays from the point of d/i formation in a pure D2, for times 
t < 20 ps and T = 30 and 300 K. 



and 300 K and <p = 0.01. The mean value of i?diff is smaller (0.80 mm at 300 K and 0.65 mm 
at 30 K, for the time interval 0-5 fis) than in the pure H 2 case since the elastic dfi(F — + d 
and dfi(F — |) + d cross sections are larger than the elastic pfi(F = 0) + p cross section 0. 
Moreover, the hyperfine splitting for dfi is AE^ = 0.0495 eV. As a result, there is practically 
no spin-flip acceleration of dfi atoms at 300 K. A relatively weaker acceleration, compared 
to the pn + H 2 case, takes place at 30 K. 

The Monte Carlo results can be compared with a simple analytical estimation. The 
kinetic theory of gases gives the following mean diffusion radius -Rdifr as a function of time: 

R 2 dm = 6Dt, (2) 

in which D denotes the diffusion coefficient. It is assumed that the atom survives until the 
time t. Using the standard definitions from the kinetic theory of gases: 

where v is the mean atom velocity, L is the mean free path, a stands for the total cross 
section and iV is the number density of atoms, one has 

— 2 i— Vt 

A»=V2— . (4) 
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The factor y/2 is valid for a simple model of the hard sphere collisions. However, the muonic 
atom scattering from hydrogenic molecules is strongly anisotropic. Therefore, we use the 
following approximation: 



where cft ra n is the transport cross section ((TJ averaged over the thermal motion of the muonic 
atoms and of the target molecules. Taking into account the muon lifetime tq = 2.2 /is, we 
obtain the following estimation of the mean diffusion radius: 



For T = 300 K, ip = 0.01 and a pure H 2 target we have a tran = 20.8 x 10~ 20 cm 2 , which 
gives i?diff ~ 11 mm. The analogous estimation for T = 30 K, using a tran = 161 x 10~ 20 cm 2 , 
leads to Rdm ~ 0.23 mm. These analytical values are in good agreement with the Monte 
Carlo results calculated assuming the thermal initial distribution of pfi energies and zero 
population of the F — 1 state (see the third column in tabled]). In the real case, the diffusion 
radius is larger owing to the epithermal diffusion. Let us note that it is very important to use 
the molecular differential cross sections for a correct Monte Carlo simulation of the thermal 
part of the diffusion. The diffusion radius in H 2 occurs to be about two times smaller than 
in the case when the corresponding nuclear cross sections are used. 

In conclusion, it has been shown that the optimal conditions for studies of the fi~ capture 
on the proton inside the ground-state pfi(F = 0) atom are achieved at the target density 
~ 0.01 when the concentration of the deuterium is depleted to the level of 10~ 7 -10~ 6 . The 
mean diffusion radius of the muonic atoms at these conditions is on the order of 1 mm. It can 
be significantly lowered when the target temperature is decreased from 300 K to 30 K. This 
effect is, however, limited as a fraction of the pfi(F = 0) atoms is epithermal both due to the 
initial high-energy component and to the deexcitation of the F = 1 states. The simulations 
of the muon capture experiments require using the differential cross sections for the muonic 
atom scattering from hydrogenic molecules. This is caused by strong molecular-binding and 
electron-screening effects at the collision energies below a feweV, where the main stage of 
the diffusion process takes place. 
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